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APPLICATION TO THE BALANCED HOMODYNE DETECTION 
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Abstract. We consider the moment operators of the observable (i.e. a semispectral measure 
or POM) associated with the balanced homodyne detection statistics, with paying attention 
to the correct domains of these unbounded operators. We show that the high amplitude limit, 
when performed on the moment operators, actually determines uniquely the entire statistics 
of a rotated quadrature amplitude of the signal field, thereby verifying the usual assumption 
that the homodyne detection achieves a measurement of that observable. We also consider, in 
a general setting, the possibility of constructing a measurement of a single quantum observable 
from a sequence of observables by taking the limit on the level of moment operators of these 
observables. In this context, we show that under some natural conditions (each of which is 
satisfied by the homodyne detector example) , the existence of the moment limits ensures that 
the underlying probability measures converge weakly to the probability measure of the limiting 
observable. The moment approach naturally requires that the observables be determined by 
their moment operator sequences (which does not automatically happen), and it turns out, in 
particular, that this is the case for the balanced homodyne detector. 



1. Introduction 

The balanced homodyne detection has long been an important tool in quantum optics, be- 
cause it is assumed to provide a way to measure the rotated quadrature amplitudes -^( e ~ Wa + 

e a*) of a single mode electromagnetic field. (See [301 EH EH El ttU E] for theoretical con- 
siderations.) An important application of this technique is the quantum state estimation, i.e. 
quantum tomography (See [23j for a collection of articles concerning this topic.) The signifi- 
cance of the homodyne measurement in that context is due to the fact that the totality of the 
rotated quadratures constitute an informationally complete set of observables, so that a quan- 
tum state is determined by their combined measurement statistics. The balanced homodyne 
detection also plays a role in e.g. continuous variable quantum teleportation [9]. 

The measurement scheme is the following: A signal field is mixed with an auxiliary field 
by means of a 50-50 beam splitter (possibly followed by a phase shifter), and the difference of 
photon numbers on the output ports of the splitter is detected. The auxiliary field is taken to 
be an oscillator in a coherent state, operating with the same frequency as the signal field. When 
the strength of the auxiliary field is high, the above mentioned photon difference statistics are 
considered to resemble those of a rotated field quadrature, the rotation angle being identified 
as the fixed phase difference between the input signal and auxiliary fields. 

The simplest description of this phenomenon is the following: Assuming that the beam 
splitter is lossless and causes no phase shift between the input modes, the process is described 
by a unitary operator U (see [21j), transforming the field annihilation operators a (signal) and 
b (auxiliary) into a = U*aU = -^(a — b), and b = U*bU = -^(a + 6). With respect to the initial 

state of the two-mode field (Heisenberg picture) , the photon number difference observable is 
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then b*b—a*a = ab* + a*b, provided that the photon detectors are ideal. (Of course, these formal 
operator relations are made precise by restricting the operators to a suitable dense subspace of 
the tensor product Hilbert space associated with the two mode system.) 

According to the heuristic explanation, the strong auxiliary field can be treated "classically", 
by replacing b with the complex field amplitude j3 = re t9 in the above operator, with 6 identified 
as the phase relative to the input signal field. Consequently, by suitably scaling the resulting 
operator with a factor proportional to the intensity of the auxiliary field, one recovers the 
rotated quadrature -^(e~ t9 a + e ld a*). 

When the auxiliary field is treated according to quantum mechanics, the high amplitude limit 
must be considered more carefully. This kind of treatment was already given in |29|, where the 
balanced homodyne detection scheme was first proposed. In that paper, however, as well as in 
[3 1 J , only the first and second moments of the photon difference statistics were calculated before 
taking the limit. As is well known, these do not necessarily determine the entire statistics. The 
limit procedure, as well as the photon difference statistics for finite amplitude, were examined 
more carefully in [8] and by considering the asymptotic behavior of the characteristic 
functions of the associated probability measures using formal operator series expansions. In 
[28j, the scheme was described in terms of positive operator measures, and the characteristic 
functions for the probability measures corresponding to the coherent input states of the signal 
field were calculated. The high amplitude limit for these functions was then shown to be equal to 
the characteristic functions for the corresponding statistics of the rotated quadratures, and the 
Levy-Cramer theorem was applied to prove that the associated probability measures actually 
converge weakly, in the sense defined in [H p. 11]. (The same theorem had been previously 
applied to the unbalanced, homodyne detector, with a similar conclusion [I7j-) However, in the 
calculations leading to the above mentioned characteristic functions, some of the problematic 
matematical questions (like term-by-term integration of a series) were not explicitly discussed, 
making it difficult to verify whether the treatment was entirely rigorous. 

In this paper, we consider the problem of the high amplitude limit in view of the moment 
operators L(x k , E) of the observables (i.e. semispectral measures) E arising from the homodyne 
detection statistics. (The moment operator L(x k ,E) is basically the "weak" operator integral 
of the power x k over E; the precise definition of this unbounded operator will be given in 
the next section.) The moment approach was used also in [3], where the authors studied the 
moment operators of the observable obtained when taking into account the imperfectness of the 
photon detectors. They viewed the moment operator sequence as the collection of "operational 
observables" describing the effects caused by the measuring arrangement, as opposed to the 
"intrinsic" observable, the rotated quadrature, which is the observable intended to be measured. 
However, they dismiss the limit problem by applying the heuristic classical approximation of 
the auxiliary field to the characteristic functions, so that their results do not contribute to our 
problem. 

We begin our approach in a more general setting: We take a sequence (_E n ) ne N of measure- 
ments (that is, observables), with the property that for each k G N and (a suitable) vector 
state (p, the numerical sequence (M* ) n6 N converges, where M* £ R is the kth moment of 
the outcome statistics obtained by measuring E n in the state ip. (Note that these moments 
are the averages of the moment operators of the corresponding observables.) Then we consider 
the question under what conditions there exists a unique limiting observable E, such that the 
moments M k of the measurement outcome statistics of E in the state p would be obtained 
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as = lim n ^oo M* . Also, we study the implications of this convergence on the level of 
the effect operators E n (B), B G <B(R), as well as the probability measures (i.e. measurement 
statistics) associated with the observables. The considerations are based on the elements of 
classical probability theory: the method of moments as well as the theory of weak convergence 
and compactness of probability measures (see e.g. [4] and [5]). 

Section [3l SI and the two appendices are devoted to the above considerations. Having investi- 
gated the general questions, we show that the moments associated with the balanced homodyne 
detection fit into this scheme: In section [5], we give a simple expression for the observable mea- 
sured by the homodyne detector as a Naimark dilation, which we use in section [6] to determine 
the essential properties of the moment operators of that observable. The purpose of section [7] is 
to demonstrate an easy (and mathematically unproblematic) way of calculating the character- 
istic functions for the probability measures corresponding to the coherent states, and thereby 
verify the similar result in [28J mentioned above. We do not, however, need these characteristic 
functions in our moment approach. 

In the final section, we gather the results together to provide a conclusion concerning the 
high amplitude limit. 

2. A FEW GENERAL NOTATIONS AND DEFINITIONS 

Let TC be a complex Hilbert space and let L(TC) stand for the set of bounded operators on 
TC. For any (not necessarily bounded) operator A in TC, we denote by D(A) the domain of 
definition of A. If A is a selfadjoint operator, we let E A : £>(R) — > L(TC) be the unique spectral 
measure of A, where £>(R) is the Borel a-algebra of R. For any two vectors (p,ip G TC, we let 
I^X^I denote the operator r\ \— > (tp\r])(p. 

If Q is a set and A a a-algebra of subsets of Q, a weakly (or, equivalently, strongly) a-additive 
set function E : A — > L(7i) with E(B) > for all B G i, and E(Q) = I, is a semispectral 
measure (also called a normalized POM). For each ip,cp G TC, the scalar valued set function 
£>(R) 3 B i — ► (if)\E(B)ip) G C is then a complex measure, and we will denote it by E^ jip . For 
each positive operator T of trace one, we let E T : A — > [0, 1] be defined by E T (B) = Tt[TE(B)]. 
Then Et is a probability measure. For a unit vector ip G TC, we then have E\ v )^\ = E V}ip , and 
we write simply E 9 to denote this probability measure. 

The operator integral [20] is the basic tool in our approach. If E : A — > L(TC) is a semispectral 
measure, and / :Q->Ca measurable function, there exists a unique (not necessarily densely 
defined) linear operator L(f, E) on TC, such that 

ty\L(f, E)<p) = j f dE^, ^eTC lV e D(f, E), 

where the domain D(f,E) is the set of those vectors ip G TC for which the function / is 
integrable with respect to the complex measure E^ }ip for all e TC. We will mostly use the 
square-integrability domain 

D(f,E) = {<peTC\ f |if d^<oo}, 

which is a (possibly proper) subspace of D(f, E) [20J. For convenience, we will use the symbol 
L(f, E) to denote the restriction of L(f, E) to D(f, E). 

When the Hilbert space TC is associated with a quantum system, any state of the system is 
represented by a positive operator T of trace one, and an observable of the system is represented 
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by a semispectral measure E : A — > £(7Y) (with the measurable space (Q,A) chosen suitably). 
The associated probability measure Et then describes the measurement statistics of E in the 
state T. In this paper, we consider only the case (Q,A) = (M, B(M)). 

The kth moment operator of a semispectral measure E : £>(M) — > L(H) is the operator 
integral L(x k ,E) (where x k is the function x \— > x fc ). Each moment L(x h ,E) is a symmetric 
operator. 

A semispectral measure i£ : B(R) — > £(7i) is a spectral measure, if E(B) is a projection for 
all £> G A Then all the moment operators L(x k ,E) are selfadjoint and densely defined, with 
D(x k ,E) = D(x k ,E) and L(x k ,E) = L(x,E) k for each k G N. Indeed, £ is then the unique 
spectral measure of the selfadjoint operator L(x,E). 

3. Convergence of observables in terms of their moments 

In the final section of the paper, we will see that the convergence taking place in the homodyne 
detection scheme can be specified by using the moment operators of the associated semispectral 
measures. The purpose of the present section is to precisely formulate the aforementioned 
moment convergence in the general setting, as well as to connect it to the weak convergence of 
the probability measures associated with the semispectral measures in question. Throughout 
the section, Tt is a complex separable Hilbert space. (Separability is needed in the main results). 

Definition 1. Let E n , n G N, and E be semispectral measures B{R) — > L(TC). If T> C 

(^n,kenD(x k , E n )) D D(x k ,E) is a dense subspace, such that 



lim {ip\L(x k , E n )y) = {ip\L(x k , E)<p), k e N, ip e H, <p e V 

n— >oo 



(in particular, each limit in the left hand side exists), then we say that E is a moment limit for 
(E n ) neN on V. 

The following observation will be useful. 

Proposition 1. Let V C TC be a dense subspace, and let E n and E be semispectral measures 
£>(R) — > L(H). Then E is a moment limit for (_E n ) nG N on V, if and only if for each unit vector 
ip G V, the probability measures E™ and E v have all moments, and 

(1) lim [x k dE"(x) = [ x k dE v (x), k G N, cp G V, \\<p\\ = 1. 
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Proof. If E is a moment limit for (E n ) nl ^ on a dense subspace V, then V C D(x 2k ,E) C 
D(x k ,E) for any fceN, and V C D(x 2k ,E n ) C D(x k ,E n ) for n, k G N, and (pp) clearly holds. 
Assume then that there is a dense subspace V, such that each x k is E™~ and E^-integrable 

whenever n G N, y? G V, \\<p\\ = 1, and that (DD) holds. Then V C D(x k ,E) C D(x k ,E), and 
the same is true for any i?* 1 in place of E. Since P is a subspace, the condition ([]]) implies that 

(2) lim (^\L(x k , E n ) V ) = (^\L(x k , E)<p). 

holds for any k G N, ip, ip G P by polarization. Let ip E V and G N. Since each £' n is a 
semispectral measure, we have \\L(x k , E n )(p\\ 2 < J x 2k dE™(x) (see e.g the proof of [2QJ Lemma 
A. 2]), and hence (pp) implies that the sequence (L(x k , E n )ip) n( z^ is bounded. It follows that (j2j) 
holds for all k G N, V> G H, ip G X> (see e.g [21 Theorem 2, p. 47]), so that £" is a moment limit 
for (E n ) nen (onD). □ 
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Remark 1. A sequence (E n ) of semispectral measures can have various moment limits E on 
a subspace V. For a trivial example, take a constant sequence E n = /if, with \i a probability 
measure which is not determined by its moments, i.e. there exists another probability measure 
with the same moments. 

As the preceding remark demonstrates, the uniqueness of the moment limit of a sequence 
of observables is connected to the limiting probability measures being determined by their 
moment sequences. Recall that a positive measure /i : £>(IR) — ► [0, oo) is called determinate, if 
each moment := J x k dfj,(x), k — 0, 1, 2 . . ., exists (and is finite) and the moment sequence 
(Mfc) determines the measure /x, i.e there is no other measure v : £>(IR) — ► [0, oo) such that 
jf x k du{x) = M k for all k = 0, 1, 2, . . .. 

Definition 2. Let E : i3(R) — > L(W) be a semispectral measure, and let V C H be a set. If 
the positive measure E v is determinate for each <p E D, then is V -determinate. 

Remark 2. It is easy to see that for a dense subspace V, a D-determinate semispectral measure 
E is determinate, i.e. the only semispectral measure having the same moment operators as 
E is E itself. In fact, assume that E is P-determinate with V a dense subspace, and let 
E' : B(M.) — > L(7i) be a semispectral measure such that L(x k ,E) = L(x k ,E') for all k E N. 
Since E v is determinate for each p E V, we have, in particular, that J x k dE v (x) < oo for all 
p E V. Hence, X? C D(x k ,E) C D(x k ,E) = D(x k ,E') for all G N. Now £^ = E' v for all 
ip E V, since ^ is determinate and L(x k , E)<p = L(x k , E')ip for all ip E V. Using polarization 
and the density of T>, we get E = E', proving that E is determinate. 

If V is not dense, a D-determinate semispectral measure need not be determinate. For 
an example, let P E L(H) be a projection, \i a determinate probability measure on R (thus 
having all moments), and define E v : B{R) -> L(H) by E V {B) = fi(B)P + v{B){I - P) 
for each probability measure v with J \x\dv(x) = oo. Now each such E v is a P-determinate 
semispectral measure, but D(x k ,E u ) = D(x k ,E u ) = P(H), and L(x k ,E u ) is the operator 
P(TC) 3 <p i— > (/ x k dfi(x)){p E H, so that the moment sequence of i?" does not depend on z/. 

Remark 3. Although the concept of D-determinate semispectral measure may seem some- 
what artificial compared to the natural definition of determinate semispectral measure, many 
existing physically relevant examples of the latter actually fall also into the former category. In 
particular, it was proved in [12] that the polar margins of a certain physically relevant phase 
space semispectral measures are determinate, but from the proof we see that these semispec- 
tral measures are actually D-determinate, where V C L 2 (IR) is the linear span of the Hermite 
functions. In [13] it was shown that the cartesian margins of the same phase space semispec- 
tral measures are C^°(M)-determinate, where C^°(R) is the space of complex valued compactly 
supported infinitely differentiable functions. 

The following result emphasizes the relevance of the concept of a "D-determinate semispectral 
measure in our considerations. 

Proposition 2. Let E n , n E N, and E be semispectral measures B(M.) — > L(7i), and V C H 
be a dense subspace, such that E is a moment limit for (E n ) n& ^ on V. If E is V- determinate, 
then E is the only moment limit for (E n ) ne ^ on T>. 

Proof. According to proposition [U any other moment limit E' on V satisfies J x k dE v = 
J x k dE'p for all <p E V. Since E is P-determinate, it follows that E v = E' for all <p E T>, 
which implies E = E' by polarization and the density of V. □ 
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It is often difficult to determine, whether a given semispectral measure is P-determinate for 
a given subspace V. In some applications, however, the relevant positive measures possess the 
following stronger property [32]: A measure \x : B(M.) — > [0, oo) is exponentially bounded, if 
J e a \ x \ d/j,(x) < oo for some a > 0. As mentioned in [12], such a measure is determinate by e.g. 
[IH Theorem II.5.2]. 

For a semispectral measure E : <B(R) — > L(TC), we let denote the set of vectors if G 7^ for 
which ^ is exponentially bounded. Later we will need the following observation. 

Lemma 1. For any semispectral measure E : £>(R) — > L(7i), the set Ee is a subspace ofTi, 
and E is Ee- determinate. 

Proof. If ip, if) G £g, and Ci, c 2 G C, we have 

E cw+C2i) {B) = \\E(B)Hci V + c 2 ^)\\ 2 < (\ Cl \\\E(B)^\\ + \c 2 \\\E(B)h\\) 2 
< 2\c 1 \ 2 E ip {B) + 2\c 2 \ 2 E lP {B) 

for all B G £>(R). This clearly implies that also E Cltp+C2 ^ is exponentially bounded, showing 
that Ee is a subspace. The last statement follows from the previously mentioned fact that an 
exponentially bounded measure is determinate. □ 

Having investigated the moment convergence in detail, we now want to describe this conver- 
gence at the level of the probability distributions themselves. The motivation for this comes 
partly from the fact that Vogel [28] claims that the relevant convergence concept in the homo- 
dyne detector problem is the weak convergence of all the probability measures associated with 
the measured observables. In the definition below, we formulate this convergence in terms of 
the observables themselves. This formulation has been used e.g. in [3, Theorem 1] (where it 
appeared without any specific name or physical context). One should recall the definition of the 
weak convergence of probability measures [Hp. 11]: a sequence (/i n ) of probability measures on 
£>(R) converges weakly to a probability measure /i : B(M) — > [0, 1] if lim n _^oo J f dfi n = f f dfj, 
for all bounded continuous functions / : R — > R. 

For each B C R, we let OB denote the boundary of B, i.e. dB — B nl\B. This notation 
is needed in a characterization of the weak convergence, which says that a sequence (fin) of 
probability measures £>(R) — > [0, 1] converges weakly to a probability measure fi : £>(R) — > [0, 1], 
if and only if \im n fi n (B) = fx(B) for all B G B(R) with /x(dB) = (see e.g. [H Theorem 2.1]). 
One should also recall that a sequence of probability measures converges weakly, if and only 
if the associated sequence of distribution functions converges pointwise at all points where the 
limiting distribution function is continuous [H p. 335-336]. 

Definition 3. Let H be a Hilbert space, and let E n : £>(R) — > L(H) be a semispectral measure 
for each n G N. We say that the sequence (E n ) converges to a semispectral measure E : £>(R) — > 
L(7i) weakly in the sense of probabilities, if 

lim E n (B) = E(B) 

in the weak operator topology, for all B G £>(R) such that E(dB) = 0. 

The following proposition characterizes this convergence. The proof of this result is given in 
Appendix A (in the more general context of a metric space). 

Proposition 3. Let E n ,E : £>(R) — > L(H), n G N, be semispectral measures. Then the 
following conditions are equivalent. 
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(i) (E n ) converges to E weakly in the sense of probabilities; 

(ii) for each positive operator T of trace one, the sequence (Ej.) of probability measures 
converges weakly to E T ; 

(iii) there exists a dense subspace V C ri, such that the sequence (£7") of probability measures 
converges weakly to E v for any unit vector ip G V; 

(iv) lim^oo L(f, E n ) = L(f, E) in the weak operator topology for each bounded continuous 
function f : R — ► R. 

The main result of this section, Proposition [5] below, gives a connection between the moment 
limit and convergence in the sense of probabilities, which is analogous to the corresponding 
result for probability measures [H2 Theorem, p. 540]. The proof of part (a) of the proposition 
is based on the following result, which uses the concept of relative compactness as it appears 
in probability theory: a family V of probability measures v : £>(R) — > [0, 1] is called relatively 
compact, if every sequence of elements of V contains a weakly convergent subsequence (see [H 
p. 35]). The proof of Proposition 0] is given in Appendix B. 

Proposition 4. Let T> C H. be a dense subspace, and let M. be a collection of semispectral 
measures E : B(R) — ► L(7i). Suppose that the set {E^ \ E G Ai} of probability measures is 
relatively compact for each unit vector tp G T>. Then every sequence of elements of M. contains 
a subsequence which converges weakly in the sense of probabilities. 

Proposition 5. Let E n : <6(R) — > L(Tt) be a semispectral measure for each n G N. Assume 
that there is a dense subspace T> C n n ^ m( z^D{x m , E n ) , such that the limit 

lim ^\L(x m ,E n )tp) 

exists in R for each m G N and <p G T>. 

(a) There exists a semispectral measure E : B(R) — > L(7i), which is a moment limit for 
(E n ) nm onV. 

(b) Suppose that E is T> -determinate. Then E is the only moment limit for (E n ) n ^ on T>, 
and the sequence (E n ) n( z^ converges to E weakly in the sense of probabilities. 

Proof. By assumption, 

sup / x 2 dE™(x) < oo 

nGN J 

for all (p G V. This implies that {E™ \ n G N} is relatively compact for each unit vector <p G T>. 
(See e.g. the first paragraph in the proof of [H Theorem 30.2, p. 408], and Prohorov's theorem 
pjj Theorem 6.1, p. 37].) Hence, we can apply Proposition [4] to get a subsequence (E nk )k£N 
converging to a semispectral measure E : £>(R) — > L(TL) weakly in the sense of probabilities. 
Fix a unit vector if G V. By Proposition [3] the sequence (E™ k )km of probability measures 
converges to E v weakly. Since, in addition, 

\ m (tp) := lim / x m dE n (x) = lim I x m dE^(x) 

n^oo J ^ k^oo J 

exists for all m G N by assumption, it follows from the proof of Theorem 30.2 of [5, p. 408] (or 
directly from [26l Theorem 4 A]) that A m (<£>) = J x m dE^(x). But by Proposition [T] this means 
that E is a moment limit for {E n ) n&i on V. Hence, (a) is proved. 

Assume then that E is P-determinate. It follows from Proposition [2] that E is the only 
moment limit for (E n ) n£N on V. Let cp G V be a unit vector. Then the numbers J x m dE v (x) 
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are the moments for only one probability measure (which is E v \ so it follows directly from [5], 
Theorem 30.2] (or [26, Theorem 6]) that (_E") ne N converges to E v weakly. By proposition [3l 
this entails that (E n ) n€ ^ converges to E weakly in the sense of probabilities. □ 

4. An "asymptotic measurement" scheme 

Suppose that we have a measurement setup which can be configured to various modes M. n , 
n = 0, 1,2, . . ., each of which is supposed to constitute a measurement, which is somehow an 
approximation of the observable we actually want to measure, the approximation becoming 
more and more accurate as n increases. This is, of course, somewhat vague, but it can be 
given a reasonably precise meaning by considering the moments of the measurement outcome 
statistics. In the following, we number the steps so that they can be referred to later. 

(1) According to quantum measurement theory, each measurement Ai n can (in principle) 
be associated with a semispectral measure E n , acting on a (common) Hilbert space 
7i. Suppose that we have a fixed set of measurement preparations to be used for the 
calibration of the measurement setup (i.e. determination of the measured observable), 
and assume that this set can be identified with the unit ball T>\ of some dense subspace 

VcH. 

(2) Then each measurement Ai n , performed with an input state ^ G Pi, produces the 
probability distribution E™, whose moments M* , k G N, can be calculated from the 
measurement data. (Notice that this can be done without needing to know the actual 
form of the semispectral measure E n ). After this point, the process obviously works 
only if we can convince ourselves that all the moments exist. 

(3) Assume then that, once a sufficient number of the measurements Ai n have been per- 
formed, one arrives at the conclusion that 

lim exists for each k G N, ip G T>\. 

Then Proposition [5] (a) says that there exists at least one observable E, such that 
the numbers obtained from the limits lim^oo coincide with the corresponding 
moments of the measurement statistics of E. Thus, in the level of moments, this scheme 
describes a valid measurement. 

(4) However, this approach cannot lead to a measurement of a unique quantum observable, 
unless the moments obtained as the above mentioned limits actually determine uniquely 
the probability distributions from which they arise. Hence, it is natural to require that 
each measure E^, ip G T>\ be determinate, i.e. E be D-determinate, in which case 
Proposition [5] (b) implies that exactly one observable E can be associated with the limit 
statistics corresponding to the set T>i of calibration states. 

(5) One additional condition could be imposed in order to make the above reasoning more 
solid. Namely, also the semispectral measures E n could be required to be P-determinate. 
In view of Proposition [5], this is not necessary, but it would nevertheless ensure that 
for each n, the moments used in the above process completely determine the 
measurement statistics of A4 n . 

Thus, if all the above conditions are satisfied, then the scheme can be considered as a mea- 
surement of the unique moment limit E of the sequence (E n ) on T>. Moreover, it is important 
to note that then (E n ) n£N converges to E weakly in the sense of probabilities (Proposition [5] 
(b)), which means (see Proposition [3]) that the sequence {EJf) of outcome probability measures 
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converges weakly to Et for any input state T (positive operator of unit trace), not just for the 
states used in the calibration. 

As we will see in the following sections, all the above requirements are satisfied for a simple 
theoretical model of the balanced homodyne detector, when the linear span of the set of coherent 
states are used in the calibration. 

Remark 4. The requirement that the observables be P-determinate is quite restrictive, and 
may be difficult to verify without knowing the explicit form for the observables. In some cases, 
the question might be resolved by using a criterion from the extensive literature concerning 
the classical moment problem (see for instance pQ or [IE]). For example, if all the limiting 

moments '■— lim^oo M* >ip satisfy the condition liminf^oo j^M^^, < oo, then the limiting 
observable E is P-determinate (see e.g. [TBI Theorem II. 5.1]). In our example of the homodyne 
detection, this condition actually holds, but we will not need to state it explicitly, as the relevant 
measures have the stronger property of exponential boundedness (see the proof of Lemma [2|). 
The fact that exponential boundedness implies the above moment condition can be seen e.g. 
from the proof of [H3, Theorem II.5.2]. 

5. The balanced homodyne detector 

Now we proceed to the description of the balanced homodyne detection scheme. A simple 
description of this detector can be given as follows: a (single mode) signal field with the Hilbert 
space TC is coupled with an auxiliary field (with the Hilbert space TC aux ) via a 50-50-beam 
splitter described by a certain unitary operator U. The auxiliary field is initialized in the 
coherent state \z). The photon number difference for the two output ports of the beam splitter, 
divided by the strength \z\ of the auxiliary field, is then detected. (We consider only ideal 
photodetectors in our scheme.) 

Assuming that TC and TC aU x are separable complex Hilbert spaces, we use the following usual 
notations. We fix orthonormal bases of the form {|n) | n G N} for both TC and TC aU x, representing 
the photon number states. We denote 

\nm) = \n, m) := \n) Cg> \m) ETC® TC aux , n,m G N. 

Let a, a* and b, b* be the creation and annihilation operators for the aforementioned bases of 
TC and TC aux , respectively, and let N = a*a and N aux = b*b, be the photon number operators 
for the two modes. The operators a, a* ,b,b* , N, N aux are considered as being defined on their 
natural domains, e.g. 

D(a) = D(a*) = {(p G TC \ ^n\(n\^)\ 2 < oo}; 

neN 



D(N) = D(a*a) = {if G TC | ^ n 2 \ (n\<p) \ 2 < oo}. 

neN 

For any z G C the coherent state \z) ETC is defined by 

00 ..n 

4N 2 



\z) = e -^ 2 ^—=\n) 



n=o v ni 



and we use the same symbols for the coherent states in TC aux . Also, we use the shorthand 
\z,w) = \z) ® \w) G TC ® TCaux- The subspace D coh := lin{|z) | z G C} is dense in TC, and 
so is the corresponding subspace D^fi in TC aux . (Here the symbol "lin" denotes the (algebraic) 
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linear span of the set in question.) The (algebraic) tensor product D^ oh := D co h <8> can 
be identified with lin{|/3, z) \ (3, z G C}, which is dense in Ti ® TC aU x- The following standard 
formulas hold: 

e itN \z) = \e u z), e itNa ^\z) = \e lt z), 
i z W\ = e -i(| z | 2 +| z '[ 2 )+^'_ 



Denote by Q and P the signal field quadrature operators ^= (a* + a) and ^= (a* — a), respec- 
tively. Here the bar stands for the closure of an operator, so that e.g. Q is the unique selfadjoint 
extension of the essentially selfadjoint symmetric operator -^(ct* + a) (See [25], Chapter IV] or 
0, Chapter 12] for details concerning the domains of these very extensively studied operators.) 
The set of states for the signal field is denoted by S(TC), i.e. S(TC) consists of those T G L(TL), 
which are positive and of trace one. 

For the signal field, define the rotated quadrature operators Qe, with 9 G [0, 2ir), via 



(3) Q e = e ldN Qe- ldN = -^{e~ ie a + e i6 a*) 

v2 

In particular, each Qe is selfadjoint on its domain D{Q e ) = e l9N D(Q) D D(a) = D(a*), and 
the restriction Qo\d{o) = '^ er%6a + e l<9 a*) is essentially selfadjoint. The ordinary quadratures 
are given by Qo = Q and Q* = P. We will need the following fact. 

Lemma 2. For all 9 G [0, 2ir), the spectral measure of Qe is D coh - determinate. 

Proof. Let 9 G [0,2-71"). By the definition ([3]) of Qe, the spectral measure E Qe is X i— > 
e l9N E®(X)e~ l6N . Let (3 G C, and consider the probability measure E9pf, which has the sim- 
ple form E Qe (X) = (f3'\E®(X)\f3'), where \(3') = \e~ l9 (3) . In the coordinate representation 
(7i ~ L 2 (R), \n) i— > h n , with h n the nth Hermite function), Q is the multiplication by the func- 
tion x I— > x, and \(3') is the vector W(q,p)ho, where (3' = ^(q + ip), and W(q,p) is the Weyl 
operator, so that (W(q,p)h )(x) = ^--^-^v+^^k^? . Hence) £^(X) = ^ J x e'^'^dx, 
which immediately implies that J e a<[X<[ dE^ < oo for any a > 0. Hence, G for each 
(3 G C It then follows from Lemma [H that is _D co /i-determinate. □ 

In quantum optics, it is typically asserted that a (lossless) 50-50 beam splitter with phase 
shift can be described by a unitary operator U^, satisfying 

(4) z) = \j= 2 ((3 - t**z\ j= 2 (e-*tf3 + z)). 

In particular, U(D co h <8> D^fi) = D co h <8> D^fi ■ In terms of annihilation operators, (TjQ) implies 

U;{a®I)\ Dlo U, = -L(a- e ^)| DLi ; 

U;{I®b)\ Dloh = -±=(e-^a + b)\ Dl 



v/2 



coh 



Using these equations, we see that the restriction of the operator I ® N aux — N ® I aux to D\ oh 
is unitarily equivalent to the corresponding restriction of {e~ % ^a®b* + e"^a* ®b), the equivalence 
being carried by U. Now both of these operators are essentially selfadjoint, and the closure of 
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the former is the photon number difference observable, which we denote by AL. Hence, 

where is the closure of the operator -^(e~ l(t> a (g> b* + e^a* ® 6). Furthermore, = I ® 

e -i<j)N a ux ^4 j (g) e #iv aTO) where ^ . = 

According to the scheme, the detection observable for the two-mode field is the photon 
difference AL, divided by the amplitude \z\ of the input \z) of the auxiliary field. For later 
convenience, we scale it with the factor so that the detection is represented by the spectral 
measure 

B(R) 3 B i-> ^ V5 '*I)~ 1JV -( J B) G L{H®H aux ) 

of the selfadjoint operator (v^D^AL. If T is the input state of the signal field, the two-mode 
field state after the beam splitter is U^T ® \z)(z\Ul. Hence, the detection statistics are given 
by the probability measures 

B^p z /(B) :=Tr[T^\z)(z\U*E^ z ^ lN -(B)U], TeS(H), z G C. 

Because of the unitary equivalences 4=[/£AL£^ = A^ = I <g> e~ i<t>Naux AI ® e i<t>Naux , we get 

(5) p z /(B) = Tr[T ® \z)(z\EW~ 1a *(B)] = Tr[T ® |e^)(e^|S w_U (S)], 

reflecting the fact that the phase shift in the beam splitter can be realized by shifting the phase 
of the auxiliary field. Accordingly, we put = in the subsequent discussions, and write 
instead z = re ld , so that 9 represents the phase difference between the two input modes. The 
beam splitter is now U := U , with 

(6) U\P,z) = \j 5 (J3-z),j-{p + z)). 

Since the signal field is considered as the input for the actual homodyne detector with the 
fixed auxiliary state \z) = \re td ), the observable E z : B(M) — > L(H) being measured by the 
detector is now uniquely determined, as a semispectral measure, by the relation 

(7) Tr[TE z (B)} = Tr[T g> \z)(z\E^ lA (B)}, B G B{R), T G S{H). 

(The reader may wish to consult [TQ] for background information on quantum measurement 
schemes and semispectral measures; the balanced homodyne detection is briefly discussed on 
pp. 194-195.) 

Let V z : TC — > TC ® Ti aux be the linear isometry <p \— > <p (g) 1^). In the formulation of measure- 
ment dilations, the relation says that (TC ® Tt aux , E^ A , V z ) is a Naimark dilation of the 
semispectral measure E z , i.e. 

(8) E Z {B) = V z *E lzl ~ lA {B)V z , B G B(R). 

(For a brief exposition of measurement dilations, see [22l Section 2].) 

The idea in the homodyne detector is that in the high amplitude limit, i.e. the limit \z\ — ► oo, 
the measurement statistics of the observables E z begins to resemble the corresponding statistics 
of the rotated quadrature Qg, where 6 is the fixed phase of the complex number z. As announced 
above, we approach this problem by considering the moments of these semispectral measures. 

11 



6. The moment operators of E z 

In this section, we determine the relevant properties of the moment operators L(x k ,E z ). 
Most of the discussion concerning the relevance of the results is postponed to the concluding 
section [3 For simplicity, we let L(x k ,E z ) denote the restriction of L(x k ,E z ) to D(x k ,E z ). 

Proposition 6. Let z G C and fc£N. 

(a) L(x k ,E z ) = \z\~ k V*A k V z . 

(b) D(a k ) C D(x k ,E z ). In particular, L(x k ,E z ) is densely defined. 

Proof. Since E^ lA is a spectral measure, we have L{x k ,E^ lA ) = \z\~ k A k , so (a) follows 
immediately from Theorem in [211 Section III A]. If (p G D(a k ), then 

V g <p = <p® \z) G D(a k ) ® D(b k ) C D((A\ D{a) ® D{b) ) k ) C D(A fc ) = P0(a^ E z ) 

by (a), so </? G D(x k ,E z ). This proves (b). □ 

The following result gives the first two moments explicitly. 

Proposition 7. Let z = re %e , with r > 0, 9 G [0, 2n). 

(a) L(x,E z ) D QelD(a)- 

(b) L(a; 2 ,^)D(g e | D(a) ) 2 + |r- 2 iV. 

Proof. Let <£> G -D(a), so that </? G D(x,E z ) by Proposition [6] (b). Let ^ G W. Part (a) of the 
same proposition now gives 

(?/>|L(a;,£ 2 V) = -^\z\-\ij® \z)\(a®b* + a* ®b)y® \z)) = (^(e^a + e i0 a*)(p). 

This proves (a). To prove (b), note first that since D(N) = D(a 2 ) = D(a*a) = D(aa*) = 
D((a*) 2 ), we have D(N) C D{{Q e \ D ^) 2 )^ and hence the domain of the operator (Qe\D(a)) 2 + 
-r~ 2 N is D(N). Now let G D(a 2 ) and E H. Using Proposition [6] again (and noticing the 
fact that (A| D(a)0l)(b )) 2 C A 2 ), we get 

{ip\L{x 2 , E z )ip) = \\z\~ 2 (il)® \z)\{a®b* + a* ®b) 2 y® \z)) 

= \\z\- 2 (i)\(-z 2 a 2 + \z\ 2 aa* + (\z\ 2 + l)a*a + z 2 (a*) 2 )^) 
= \{i)\(e- 2l0 a 2 + e 2ie (a*) 2 + aa* + a*a + r" 2 a*a)^) 
= ^|([^(e-^a + e l V)] 2 + ir-Va)^). 

□ 

Corollary 1. The intrinsic noise operator Af(E z ) := L(x 2 ,E z ) — L(x,E z ) 2 of the observable 
E z is the selfadjoint operator \r^N . In particular, D(L(x 2 ,E z )) fl D(L(x, E z ) 2 ) = D(a 2 ). 

Proof. Since D(N) C D((Q e \ D{a) ) 2 ), we have D(N) C D(L(x, E z ) 2 ) by part (a) of the preceding 
proposition. In addition, D(a 2 ) C D(x 2 ,E z ) by part (b) of that proposition. Hence, D(N) C 
D(L(x 2 ,E z ) - L(x,E z ) 2 ), so that C L(x 2 ,E z ) - L(x,E z ) 2 by the proposition. But 

L(x 2 ,E z ) — L(x,E z ) 2 is symmetric (even positive), and \r^N is selfadjoint, so that \rbN = 
L(x 2 ,E z ) - L(x,E z ) 2 . □ 
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Remark 5. Proposition [7] (a) reflects the fact that no E z is a spectral measure (for any spectral 
measure is well known to have zero intrinsic noise). A perhaps somewhat less well-known fact is 
that a semispectral measure with zero intrinsic noise, and a selfadjoint first moment operator, is 
necessarily a spectral measure, and thereby completely determined, as a semispectral measure, 
by its first and second moment operators (see pEH Theorem 5]). By the above corollary, we 
have J\f(E z ){p — ► as \z\ — > for any ip G D(H{E Z )) = D(N), so one might expect that the 
first two moments of the detection statistics indeed determine the observable measured in the 
high amplitude limit. However, due to the problems concerning the domains of the moment 
operators, we were not able to make this reasoning rigorous, and therefore chose instead to 
consider entire moment sequences in the "weak sense", as will be shown in the last section. 

Remark 6. In view of the tomographic application, there is also another aspect in the above 
mentioned fact that it is not enough to consider only the first moment operator. Namely, 
one cannot reconstruct a state by knowing only the averages of all the quadratures in that 
state. As a simple example, consider any two different vector states \<p)(ip\, for which 

<f,ip G D(Q) H D(P), (<f\Q<f) = (iplQip), and ((p\P(p) = (ip\Pip). [For example, one can take, 
in the coordinate representation, ip(x) = f(x)e l ^ x \ ijj(x) = f(x)e~ % ^~ x \ where / is a positive 
even function, and <fi : R — > [0, 2n) such that x h- > 6{x) + 9(—x) is not a constant mod 2ir.] 
Since D(Q) D D(P) = D(a), it follows that (<p\(up) = {^aip) and (<p\a*ip) = (ij;\a*i/;), and, 
consequently (ip\Qgip) = (ip\Qoip) for all 9 G [0, 2tt). Hence, the averages do not distinguish 
between these states. 

For the general moment operator L(x k , E z ), the following result is sufficient for our purposes. 
Proposition 8. Let k G N, and z = re ie , with (r, 9) G [1, oo) x [0, 2tt). Then 

L(x\ E z )\ D(ak) = (Q e \D(a)) k + ^C k (r, 9), 

where 

C k (r,9)= c k nm (r,9)(aTa m , 

n+m<k 

and each function c k im : [1, oo) x [0, 2ir) — > C is bounded. 

Proof. Let (r,9) G [l,oo) x [0,2tt), z = re ie , and let ip G H, (p G D(a k ), so that <p ® \z) G 
D(a k ) ® D(b k ) C D(A k ). By using again Proposition El we get G D(x k 1 E z ), and 

(?P\L(x k ,E reie )p) = (i;®\z)\{{V2r)- x A) k ip®\z)) = 2^r~ k {^® \re ie )\{a®b* + a* ®b) k p®\re ie )) . 

The expression (a®b*+a*®b) k is the sum of all products of the form A\Ai ■ ■ ■ A k ®BiB 2 ■ ■ • B k , 
where for each i = 1, . . . , k, either A4 = a and E>i = b* , or Aj = a* and E>i = b. We can write 
a particular product as A1A2 ■ ■ ■ A k <8> BiB 2 ■ ■ ■ B k = p(a, a*) ® p(b*, b), where the symbol p(-, •) 
represents the rule telling which Ai'.s coincide with the operator in the first argument, and 
which with the one in the second argument. Let U k denote the set of all such rules p, and let 
C n fc be the set of those p for which the operator in the first argument appears exactly n 
times in the resulting product. Now Ii k = U,^ =o n^, and we have 

k 

(a®b* +a* ®b) k = ^2^2 P( a > a *) ® P( b *> h )- 
™=o pG n£ 

13 



Therefore, 

k 

(ij\L(x k ,E^ 9 )v) =2-*Y; Y,^\P^ a *)v){z\r- k v(b\b)\z). 

n=o pen* 

For p G n„, the operator r~ h p(b* ,b), when brought to the normal order, is of the form 
r~ k p(b*, b) = r - k (p*y i l) k - n + r~ k Rp(b*, b), where R* is a polynomial of b* and 6 of order k — 2, 
with the coefficients depending on p. Note that in each term of R* the number of occurrences 
of b* minus the number of occurrences of b is the same, namely n — (k — n). (ic is given by 
the Wick rule.) Hence, 

(z\r- k p(b*,b)\z) = (e- ie ) n (e ie ) k - n (l + r- 2 d k p (r)), 

where each function dz : [1, oo) — > M is bounded. (Of course, we could have taken any interval 
[5, oo) with 5 > 0.) For p G IT*, we have (e" ie ) n (e i9 ) fe ~Xa, a*) = p(e~ id a, e w a*), so that 

n=o pen* 
fe 

n=o pen* 




with 

tf fc (r,0) = 2-1 £ 4(r)p(e- l Ve l V). 
pen fc 

By bringing each p(e~ ld a, e ld a*) into the normal order, and using the fact that each d k is a 
bounded function, we see that Ck{r,9) has the required form. □ 

The following result shows, in particular, that each observable E z is uniquely determined by 
its moment operator sequence (L(x h , E z ))k£N (see Remark[2]). Hence, it is possible to consider 
the moment sequences in place of the observables themselves. 

Proposition 9. For each zGC, the semispectral measure E z is D ' co ^- determinate. 

Proof. Let z G C. Now D coh C D(a n ) C D(x n ,E z ) for each n G N, so for each <p G -D coft , the 
positive measure E* has all moments. 

Consider the positive measure E z ^, for fixed (3 G C. According to (JBJ), we can write Z7|/3, z) = 
\c,d), with c = i(/3 - z), d = + z). Since A = -j=U*N_U, it follows from © that 

i?^ = i?^' 2 ''' 7V_ . This measure is supported in (v^l-zl)^ 1 ^, with density 

(v^M)- 1 ^ Yl \(c,d\n h n 2 )\ 2 , 

rti,n 2 eN, 
n\—n,2=k 
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so that 

e o|x| dEj^(x) = e a(v/5|2|rl|ni - n2| | 



[M^Udln^ 



-(|c[ 2 +|d| 2 ) e a(v / 2kl)- 1 K+"2) l£l Pi 



< e -(|c| 2 +|rf| 2 ) ^ 



nr. no! 

ni,n 2 =0 

= exp[-(|c| 3 + |d| a ) + e^'l^dcl 2 + |d| 2 )] 
= exp[(|/3| 2 + |^| 2 )(e a(v/I|2|)_1 - 1)] < 00. 
Hence, \(3) G Ee* for all j3 G C. It remains to apply Lemma Q] to complete the proof. □ 

7. A REMARK CONCERNING CHARACTERISTIC FUNCTIONS 

Before proceeding to the conclusion of the paper, we calculate the characteristic functions 
of the probability measures of the observables E z , associated with the coherent states. The 
purpose of this is to verify that the results are the same as those of Vogel [28], obtained 
by a different method. Also, we want to check that the moments calculated by using these 
characteristic functions are consistent with the moment operators we obtained in the preceding 
section. 

Define $ z (t) = L[e lt \ E z ) for each t G K. Since the function x 1— > e ttx is bounded, each $ z (t) 
is a bounded operator. Hence, by using jEJ) and the Theorem in Section III A], we have 

$(t) = V z *e u ^~ lA V z = V z *U*e lt ^ zl) ~ lN -UV z . 
By using ([6]) and putting c = — z ), d = ^75 (/^ + z), we then get 

Mt) ■■= mm 

= (cle-^kD-^W^Igi^N)- 1 ^ = e -\c\Hce-^^r\ e -\ d \2 + de^^\)- 1 d 



|2 \o\2 



exp (_(| C |2 + | d |2) + e -itm Z \)-^ c \2 + e ft(V2[z|)-i| d |i 

+ §(/3 -z)(fi- zy-^VAr 1 + + z )(p- z)e u ^ z ^ . 



exp 



When replacing z with iz (i.e. shifting the phase by |), and t with —t (accounting for the 
fact that we used / <g> N 2 — Ni <S> / for the photon difference, instead of N% <g> I — I <%> iV 2 ), and 
omitting the scale factor a/2, we can compare the above result with that of Vogel [28, p. 17], 
and we see that they are indeed the same. 

Now <fp(t) should generate the moments of the probability measure EL, via the formula 



\L(x\E z )\[3) = j x k dEfo=C k ^(0). 



We check the first two moments for consistency. Let again z = re td . Clearly, (fp(0) = 1. 
Differentiating ipp gives 

<Pp(t) =i^k|- 1 (-|c| a e^^- 1 ' + |d|V^- 1 *)^(t), 
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so that <p'p(Q) = i^\ z \ 1 (~\ c \ 2 + \d\ 2 ) = ~^\ z \ 1 ( Z P + z P)i an d hence the first moment of the 
probability measure Efo is 

-i^(0) = ^(e~ ie /5 + e ie p) = (/3|i(e- i0 a + e»a*)|/?>, 

in agreement with Proposition [7] (a). 
In addition, 

^( t ) = _i| 2 |-2 ((IclV^I^" 1 * + Id^e^W) -1 *) + (_| c | 2 e -*(^l'D- 1 * + Idpe^l*!)" 1 *) 2 ) <^(t), 
so the second moment is 

rVj(o) = |N- 2 (|/5| 2 + H 2 ) + i|^|- 2 [(^ + ^)] 2 

= |([(e- l9 /3 + e^)] 2 + l) + i|/5| 2 r- 2 
= ±(e" 2 ^ 2 + e 2 ^ 2 + (|/3| 2 + 1) + \[3\ 2 ) + ±|/3| 2 r" 2 
= (/3|±(e~ 2i V + e 2i V) 2 + ««* + a*a)\/3) + (/3||r- 2 a*a|/3> 
= + A*)] 2 + §r- 2 a*a)|/?>, 

which agrees with the equation of Proposition [7] (b). 

8. Convergence in the high amplitude limit 

We are now ready to establish a conclusion on the high amplitude limit, by connecting the 
results concerning the detector observables E z with the "asymptotic measurement" procedure 
described at the end of Section [3j Since the considerations of that section involved sequences of 
semispectral measures, we need to fix a sequence (r n ) of positive numbers converging to infinity. 
For this choice, let z n {6) = r n e l6 , where the phase 6 G [0, 2tt) is also fixed. 

Choose the set D co h of coherent states to be the subspace containing the calibration states 
mentioned in part (pQ) of the aforementioned procedure. It follows immediately from Proposition 
[8] that for each 6 e [0, 2%), the spectral measure of the rotated quadrature Qg is a moment limit 
for (E Zn ^ d ') ne ^ on D co h, so that the requirements for parts (J2J) and j3]) are satisfied. The 
determinacy requirement of part (J4|) is given by lemma so that Proposition [5] (b) implies that 
the sequence (E z "^) neN has only one moment limit on D coh , namely E Qe . 

Finally, the requirement for part (jSJ) is satisfied by proposition [9l Thus, we can conclude 
that the (above simple theoretical description of the) balanced homodyne detector fits into our 
"asymptotic measurement" scheme. In particular, the sequence (E Zn ^) n&i of observables con- 
verges to the spectral measure of the rotated quadrature Qg weakly in the sense of probabilities. 

Remark 7. By definition, the weak convergence in the sense of probabilities of (E Zn ^) neN 
to E Qe means that the sequence (E Zn ^ (B)) ne jn converges to E Qe (B) in the weak operator 
topology, provided that E Qe (dB) = 0. Since Qg is unitarily equivalent to the position operator 
on £ 2 (IR), this condition is equivalent to X(dB) = 0, where A is the Lebesgue measure of the 
real line. In particular, all intervals (finite or infinite) qualify. 

Remark 8. It should be emphasized that lim n _ >clo £' 2! "W(S) = E Qe (B) (in the weak operator 
topology) is not true for all B e Z3(R). In fact, let B = |J ngN (v / 2r n ) _1 Z, where z n {9) = r n e %e 
as before. Now B is countable, and, in particular, B 6 <6(R). It follows from (jHJ) and the 
relation A = -^U*N-U that that each semispectral measure E Zn ^ is supported on (v / 2r n ) _1 Z. 
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Hence, E Zn ^(B) = I for all n G N, so that (-E 2n ^)(5)) neN converges to / in the weak operator 
topology. But E® e (B) = O, because B is countable and E Qe is unitarily equivalent to the 
position operator on £ 2 (R), thereby having the same null sets as the Lebesgue measure. On 
the other hand, B is dense in R, so dB = R, and consequently E® e (dB) = /, which is not zero. 

The above conclusions are summarized as follows. 

(a) The observable E z measured by the balanced homodyne detector is -Dc^-determinate 
for each z G C; 

(b) The spectral measure E Qe of the rotated quadrature Qe is -Dco^-determinate for each 
#g[0,2tt). 

(c) For each 9 G [0, 2tt), and each discretization z n {6) = r n e ld , lim^oo?^ = oo, the spectral 
measure E® B is the unique moment limit of (E Zn( - 9 ^) n( z^ on D co h. 

(d) The sequence (E Zn ^) ne ^ converges to E Qe weakly in the sense of probabilities. 

(e) For each state T G S(H), we have lim^^ Tr[TE*»W(B)] = Tr[TE Qe (B)], whenever 
X(dB) = 0, where A is the Lebesque measure on R. 

Appendix A 

In this appendix, we give a proof for proposition [3l Since it does not depend on the special 
properties of the real line, we state it here in the context of a general metric space. For the rest 
of this section, we fix fl to be a metric space, with metric d, and let B(O) denote the associated 
Borel cr-algebra. As before, TC is a complex separable Hilbert space, and the general notations 
given in section [2] will be used. 

We recall the definition of the weak convergence of probability measures pfl, p. 11]: a sequence 
(fjL n ) of probability measures on 13 (fl) converges weakly to a probability measure /i : J3(Q) — > 
[0, 1] if lim^oo J f dfi n = J f dfi for all bounded continuous functions / : fl — > R. We will use 
the following characterization for the weak convergence [4J Theorem 2.1, p. 11]: sequence (//„) 
of probability measures on B(fl) converges weakly to a probability measure /i : J3(Q) — > [0, 1] if 
and only if lim^oo fi n (A) = fi(A) whenever A G B(Q) is such that fi(dA) = 0. At the level of 
semispectral measures E : B(Q) — > L(H), a natural analogue for the weak convergence is the 
following. 

Definition 4. Let E n : B(Q) — > L(H) be a semispectral measure for each n G N. We say that 
the sequence (E n ) converges to a semispectral measure E : B(Q) — > L(H) weakly in the sense 
of probabilities, if 

lim E n (A) = E(A) 

in the weak operator topology, for all A G B(Q) such that E(dA) = 0. 

The next proposition characterizes this convergence in terms of the weak convergence of 
probability measures. First we need the following lemma. 

Lemma 3. Let V be an at most countable collection of finite positive measures v : B(£l) — > 
[0, oo). Define T v = {A G B(Q) | v(dA) = for all v G V}. 

(a) JF-p is an algebra which generates the Borel a -algebra B(fl). 

(b) If (j2 n ) is a sequence of probability measures on B(Vt), such that lim^oo fi n (A) = n(A) 
for all A G T-p, where /i : B(fl) — > [0, 1] is also a probability measure, then (/z„) converges 
to ji weakly. 
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Proof. Since d(A U B) C dA U 95, and <9A = <9(fi \ A) for any i,5cfi, it is clear that .F-p is 
an algebra. The rest of the proof uses an argument similar to one appearing in the proof of [27, 
Theorem 4]: Let G C Vt be an open set. For each 5 > 0, let A s = {x G | d(x, £l\G) > 5}. 
Now Ag is open, and hence a Borel set. Since OAs C {x G O d(x,Q \ G) = 5}, the sets 
dAg are disjoint for distinct values of 8. Let v E V. Since ^ is a finite positive measure, the 
family {v(dAs) \ 5 > 0} is summable, and hence the set {5 > | //(cM,;) > 0} is at most 
countable. Since V is at most countable by assumption, also the set {5 > | Ag Jy>} has 
this property. This implies that there exists a sequence (5 n ) of positive numbers converging to 
0, with A$ n G T-p for all n EN. But clearly G = UneN^<5n7 so 2~v generates B(Q). This proves 
(a), and (b) follows from [U Theorem 2.2, p. 14]. □ 

Some parts of the proof of the following proposition can be extracted at least from [7] or 
[27] (the latter concerning only multiplicative operator measures), so one might expect that the 
result is essentially well-known. However, having not been able to find the entire proof in the 
literature, we give it here. 

Proposition 10. Let E n : B{Q) — > L(H) be a semispectral measure for each n G N, and let 
also E : B{VL) — > L(H) be a semispectral measure. Then the following conditions are equivalent. 

(i) (E n ) converges to E weakly in the sense of probabilities; 

(ii) for each positive operator T of trace one, the sequence {Elpj of probability measures 
converges weakly to Et; 

(iii) there exists a dense subspace V C 7i, such that the sequence (E™) of probability measures 
converges weakly to E v for any unit vector ip G V; 

(iv) limbec L(f, E n ) = L(f, E) in the weak operator topology for each bounded continuous 
function f : Q — > R. 

Proof. To prove that (i) implies (ii), we assume (i), so that the sequence (((p\E n (A)(p)) n( zfq of 
numbers converges to (<p\E(A)(p) for each ip G 7i and A G Q where Q = {A G B(Q) \ E(dA) = 
0}. Since E is a semispectral measure and H is separable, there exists a finite positive measure 
v : B(Q) — > [0, oo) with the same sets of measure zero as E. (One can take, for instance, 
v(A) = Tr[SE(A)], with S = J2n&N ^IVn) (Vn\, where {r] n \ n G N} is an orthonormal basis of 
TC.) Hence, using the notation of Lemma [3] we have Q = T{ v }. 

Now fix a positive operator T G L(H) with unit trace, and let A G Q. Write T = 
XlfcLi tk\(Pk)((Pk\, where the <pk are unit vectors and Yln^n = 1; the series converging in the 
trace norm. Now we have 

oo oo 

lim Ty[TE u (A)} = lim V t k (ip k \E n {A)^ k ) = V ' t k lim {<p k \E n {A)<p k ) 

n^oo n— >oo ' ' ' « n— »oo 

k=l k=l 

oo 

= ^2t k ((p k \E{A)ip k ) = Tt[TE(A)], 

k=l 

where the change of the order of the limit procedures is permissible, since the inequality 
||-E n (^L)|| < 1 implies that the series 

Y,h{Vk\E n {A)<p k ) 

k=l 
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converges uniformly for nGN. Hence, (E^(A)) converges to Et(A) for all A E Q. Now Lemma 
[3] (b) gives (ii). 

Clearly (ii) implies (iii). Now we assume (iii) and prove that (i) holds. Let V be the dense 
subspace in question, and A E B(Q) be such that E(dA) = 0. Then E^(dA) = for all unit 
vectors >p E T>, so that lim^oo E™(A) = Ep(A) for all y? G P, ||</?|| = 1 by (iii). Let EH, and 
select a sequence ((p^) converging to <p, such that p> k ED for all fc. Then 

lim (p\E n (A)p) = lim lim (<p k \E n (A)p k ) = lim lim (p k \E n (A)p k ) = ((p\E(A)(p), 

where the change of the order of the limit procedures is justified since the k limit is uniform in 
n. Hence, (i) follows by polarization. 

Finally, it is clear from [H Theorem 2.1] that (ii) implies (iv) and (iv) implies (iii), since 
e.g. (<p\L(f, E)(p) = J fdE v for any unit vector <p E H, and a bounded continuous function 
/ : Q — > R. (Notice that e.g. L(f, E) E L(H) because / is bounded.) The proof is complete. 

□ 

Remark 9. Since a sequence of probability measures cannot converge weakly to two different 
limits [H Theorem 1.3, p. 9], the above proposition shows, in particular, that a sequence 
of operator measures can converge to at most one operator measure weakly in the sense of 
probabilities. 

Appendix B 

This appendix is devoted to the proof of Proposition [H We use the same general context of 
a metric space as in Appendix A, as well as the definitions and notations given there. 

In the context of probability theory, a family V of probability measures v : B(Q) — > [0, 1] 
is called relatively compact, if every sequence of elements of V contains a weakly convergent 
subsequence (see [4, p. 35]). 

It is an old result that for a sequence of probability measures on R, the convergence of 
moments (which was a central concept in the considerations of section [3|) and weak convergence 
are related, in the case where the limiting measure is determinate (see [5], p. 405-408], the 
original paper [15] by Frechet and Shohat, and also [26]). This is connected with the fact that 
the relevant sequences of probability measures are relatively compact (see e.g. [U Theorems 
6.1 and 6.2, p. 37] and the proof of P, Theorem 30.2, p. 408]). The following result shows 
that the relative compactness of probability measures is reflected in a natural way in the level 
of operator measures. 

Proposition 11. Let T> C Ti be a dense subspace, and let Ai be a collection of semispectral 
measures E : 13 (CI) — > L(7i). Suppose that the set {E v \ E E M.} of probability measures is 
relatively compact for each unit vector ip E V. Then every sequence of elements of M. contains 
a subsequence which converges weakly in the sense of probabilities. 

Proof. Let {tp n \ n E N} be an orthonormal basis of Ti, included in V. (Since TL is separable 
and V is dense, such a basis exists.) Let V be the linear span of this basis. Then let 

/C = {||«Vn + (Pm\\~ l {i l( Pn + Vm) | U, 771 G N, / = 0, 1,2,3}. 

Now fC is countable, so we can put /C = {ipj \ j E N}, where ipi ^ ipj f° r i 3- 

Let (E n ) be a sequence of elements of M.. For each j E N and n E N, let E™ denote the 
probability measure E^, . Since ipj E V, the sequence (E") ne s (and every one of it subsequences) 
contains a weakly convergent subsequence by assumption. 
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Now we can use the classic "diagonal process" to get a subsequence (E nk ) of (E n ), and a 
sequence of probability measures (fJ>j), such that (Ej k )k<^m converges to fij weakly for all j G N. 

The rest of the proof uses somewhat similar arguments to the proof of [3 Theorem 1]. 

For each i,j G N, I = 0,1,2,3, define := fi n , with n G N (the unique index) such that 
ip n = \\i l <fi + ipj\\~ 1 (i l ipi + (fj). Then for each i,j G N, define 

1 3 

1=0 

Now each is clearly a complex measure, with the property that Uij(fl) = (<Pi\<Pj) for any 
Finally, for each ip, tj) G V, set 

(9) G(A)(ip, <f) — Y,{^Wi){<PiW>iA A l AeB(Sl)^^eV. 

i,j<=N 

Note that the above double sum is always finite, since (p, ip G V. Hence, the set function 
B(Q) 3 A i— > G(A)(ip,<p) G C is a complex measure, which satisfies 

(10) G(n)ftM = <Vb>, Vwev. 

For each A G B(Q), we now have a map G(A) : V x V — > C, which is clearly sesquilinear by its 
definition ((9j). We have to prove that it is positive, i.e. 

(11) G(A)((p, <p)>0, <p G V, A g B(n). 
To that end, let V = {/ij \ j G N}, so that 

lim (ipAET^A)^) = N (A), 3 eN,Ae T v , 

k—foc 

with T-p as in Lemma [3l This means that 

]iia(i l cp i + cp j \E nk (A)(i l cp i + (p j )) = ¥ l Vi + Vj\\ 2 Ai{A), i,j G N, 1 = 0,1,2,3, AeJ 7 ?. 

k-^oo 

It now follows from ([9]) (by using the polarization identity to the sesquilinear forms (if), ip) i— > 
(V>|£ nfc (AV)) that 

(12) lim(V>|£ nfe (AV) = G(A)(V>,^), Ae^^eV. 

Since each E Uk (A) satisfies < (ip\E nk (A)ip) for all ip G V, this clearly implies that the same is 
true for G(A), provided A G T-p. We have to prove that this is true for all A G B(Q). To that 
end, let 

B x = {A G B(Q) | < G(A)(<p, ^) for all y> G V}. 

Now jFp C B\. We show that i5i is a monotone class. Let (An) be an increasing sequence of 
sets in B\, and let A = U^ =l A n . Let ip G V. Since -B i— > G(B)(<p,<p) is a complex measure, we 
get 

G(A)(<^) = lim G(A0 >0, 

n— >oo 

so that U^ =1 A n = A G B\. Similarly, we see that P^ =1 A n G B\ for a decreasing sequence (An) 
of sets of B\. This shows that B\ is a monotone class. Since T-p C it follows from the 
monotone class theore that B\ contains the a-algebra generated by T-p. But this is B(Q) by 
Lemma [3] (a), so that ffTTjl holds. 
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Now (PUD and (PSJ) imply that < G(A)(p,y>) < |M| 2 for any ^ e V and A e By 
polarization, this implies that 

sup{|G(A)(^V)| I ^ e e V, |M| < 1, |M| < 1} < oo. 

Since V is dense, it follows that for each A e B(Q), there exists a unique operator E(A) G L(TC), 
such that G?(A)(y7,V) = (vl-^(^)^) for all ip,ip e V. Now ([II]) and (PUJ) imply that the map 
E : B{Q) — > L(7i) is a semispectral measure. 

By CE2D, hm^oo E^(A) = E V (A) for all A e T v and each unit vector (p e V. It follows from 
Lemma [3] (b) that for each unit vector <p 6 V, the sequence (E™ k )kem of probability measures 
converges weakly to the probability measure According to Proposition [101 this means that 
(£' nfe ) fcgN converges to weakly in the sense of probabilities. □ 
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